Giant infrared intensity of the Peierls mode at the neutral-ionic phase transition 
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We present exact diagonalization results on a modified Peierls-Hubbard model for the neutral-ionic 
phase transition. The ground state potential energy surface and the infrared intensity of the Peierls 
mode point to a strong, non-linear electron-phonon coupling, with effects that are dominated by 
the proximity to the electronic instability rather than by electronic correlations. The huge infrared 
intensity of the Peierls mode at the ferroelectric transition is related to the temperature dependence 
of the dielectric constant of mixed-stack organic crystals. 

PACS numbers: 78.30.Jw, 63.20.Kr, 74.25.Kc 
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The nature of electron-phonon (e-ph) coupling in corre- 
lated electron systems is an enduring topic in condensed 
matter research, currently as experimental indications of 
strong, anharmonic e-ph coupling have been recognized 
in high-T c superconductors based on CuO jjj or the novel 
MgB2 0. The topic itself is old and deeply rooted in 
the problem of structural deformation in systems close 
to electronic instabilities, where particularly high sensi- 
tivity to e-ph coupling is expected. Low-dimensional sys- 
tems have a host of electronic instabilities such as Mott, 
Peierls, and spin-Peierls transitions [|| with characteristic 
vibrational signatures [Q. All these instabilities are ac- 
companied by lattice relaxation representing the extreme 
consequence of the non-linearity of electronic responses 
to perturbation by phonons. 

The neutral-ionic transition (NIT) in a rigid donor- 
acceptor (DA) chain is the boundary between a param- 
agnet and a diamagnet ||, analogous to a Mott metal- 
insulator transition. Such electronic instability is again 
accompanied by structural relaxation. NITs were ob- 
served in organic charge-transfer (CT) crystals with 
mixed DA stacks ||, and, as recognized from the out- 
set Q, dimerization accompanies the transition, leading 
to ferroelectric (FE) behavior || for proper interstack ar- 
rangement. Recently NIT models were proposed for FE 
oxides Q . In this context the anomalous role of dimeriza- 
tion in displacing electronic charge has been underlined, 
to suggest enhanced e-ph coupling near the NIT. This 
was confirmed by subsequent calculations [ fl0| without, 
however, pursuing experimental implications. 

In this Letter, we exploit the recent definition [ll[] 
of polarization, P, in insulators with periodic boundary 
conditions (PBC) to obtain the infrared (IR) intensity of 
the dimerization mode. The IR intensity depends sensi- 
tively on structure, which we relax fully, and peaks at the 
FE transition of either correlated or uncorrelated models, 
clearly implicating the structural instability rather than 
correlations for the huge intensity. Working with a lin- 
ear and harmonic model for e-ph coupling, we recognize 



large non-linearity in the anharmonicity of the ground 
state (GS) potential energy surface (PES), and within 
the adiabatic approximation we demonstrate an intimate 
entanglement of electrons and phonons. Sizeable non- 
adiabatic effects are expected near the NIT 12 , but the 



subtle interplay between electrons and phonons we dis- 
cuss here sets the baseline for identifying non-adiabatic 
effects. 

We describe NIT in terms of a Peierls-Hubbard model 
with alternating on-site energies (A) and one electron per 
site. As in the Su-Schrieffer-Heeger (SSH) model jl3|, lin- 
ear Peierls coupling to a k = phonon that modulates 
the transfer integrals is introduced in the following adia- 
batic Hamiltonian: 



H = -£(!+ (-ITS) (cLc m , ff + h.c.) (1) 
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where N is the number of sites, c\ creates an electron 
with spin a on site i, and the last term accounts for 
the elastic energy. Accounting for alternating charges on 
molecular cores, qj=2, on D, A sites, respectively, we 
define the ionicity operator as p = J2j(lj ~ hj)/N. 

For A = Eq. (@) reduces to the SSH model || in 
the U = limit, or to a spin-Peierls model in the U » 1 
limit, and the Peierls instability leads to finite dimeriza- 
tion amplitude, S, in either case. The I phase with one 
electron per site (p ~ 1) is always dimerized jl4|, p"5| . 
With increasing A, electrons are paired on —A sites in 
the N phase and the GS is a regular chain (5 = 0). The 
equilibrium 5 decreases with A and vanishes for A larger 
than a critical value that depends on the lattice stiffness, 
as shown by Rice and Mele for U = p6[ . The correlated 
case is a Mott insulator on the I side [O , with an NIT to 
a band insulator around A ~ {7/2 in the regular chain. 
Accordingly, we vary r = A — {7/2 in the vicinity of the 
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FIG. 1: The GS energy of a 12 site ring with e d = 0.64. 



FIG. 2: The fluctuating charge contribution to dP/dS, cal- 
culated for N = 12 and U = oo. Continuous line refers to the 
regular chain. Dotted and dashed lines refer to the equilib- 
rium dimerization starting at the intersection with the con- 
tinuous line. The inset shows results for U — and N = oo. 



NIT and solve Tt(T, ed) exactly in the limit A, U — > oo for 
finite chains with up to 12 sites and PBC, as required for 
the P operator ML using the diagrammatic valence bond 
(VB) technique |lq| . The lengthy geometry optimization 
limits N, but is compensated by the rapid convergence 
of GS properties compared to excitation energies |lj ■ 

In the rigid lattice (ed = 0) p varies smoothly and the 
N-I crossover, at T c = —0.668 and p ~ 0.63, is signaled by 
the unconditional dimerization instability of the I phase 
|jr7f . A deformable lattice (finite ed) dimerizes on the N 
side Jlf|. Fig. [j] depicts the evolution with V of the GS 
PES of a 12-site ring with ed = 0.64. To focus attention 
on dimerization, the GS energy for the regular chain is 
set to zero at each T. Far in the N phase (large posi- 
tive r) the regular stack is stable, but with decreasing 
r, dimerization sets in well before r c , at p ~ 0.25. In 
the pre-transitional regime the GS PES is clearly anhar- 
monic. The two minima in the I regime correspond to 
ferroelectrics with opposite polarization. 

To calculate the IR intensity of the Peierls mode, pro- 
portional to (dP/dS) 2 , we slightly modify the Resta's 
twist operator |ll| , to account for alternating charges on 
molecular cores, as follows: 



P = ^-Im(ln(^\e l 
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where "J/ is the GS wavefunction, and M = ■ Tj(qj—hj), 
with Tj locating the j site. The Af-operator is diagonal 
in the VB basis and the calculation of P is easily imple- 
mented. M accounts for the actual geometry of the stack 
and depends on 8 according to: 



- N 2 m-^ „ N8 „ 

M = > irii H n 
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where the mean intersite distance and the electronic 
charge are set to 1. The dimensionless e-ph coupling 



constant, a — y/Ked, is estimated ~ 5-10, for typical 
values of the elastic constant, K. dP/dS can be rigor- 
ously decomposed into two contributions, one obtained 
by allowing M to vary with 8, at fixed 'F, the other one 
by varying 'F at fixed M. The first contribution describes 
the IR intensity originating from the displacement of sites 
carrying a fixed charge, p. Frozen charges account for IR 
intensity of optical modes in ionic lattices. The second, 
electronic contribution instead accounts for the charge 
flux driven by 6. It only appears in systems with mobile 
or delocalized electrons. 

For a regular chain, 8 — 0, the reflection mirror residing 
on each site ensures real Z values Ej . Z is positive in 
the N regime, negative in the I regime, and Z = locates 
the NIT of a 8 = ring of any size with PBC 0. The 8 
derivative of P in Eq. (||) then diverges at the NIT, where 
the GS is metallic [0 . In dimerized chains Z is complex 
and the divergence of dP/d8 is suppressed. Since the GS 
of Eq. (1) always dimerizes on the N side of the NIT, the 
IR intensity remains finite for the equilibrium structure. 
The IR intensity of the regular chain is then an upper 
bound for the relaxed structure that, depending on ed, 
holds on the N side of the FE transition. 

We numerically evaluate dP/d8 at the equilibrium 5 
and distinguish between contributions due to fluctuat- 
ing and frozen charges. The latter is small in the en- 
tire range of T. Large IR intensity is associated with 
fluctuating charges. Fig. || shows (8P/d8)fi for 12-site 
chains. The 8 = curve diverges at T c on both the 
N and I side. Giant IR intensity at the NIT is due to 
finite electronic fluxes generated by infinitesimal lattice 
distortions. Resta and Sorella reported similar results 
by finding the total average charge for fixed dimerization 
JlOj l : for small 8, their quasi-divergent curves are finite- 
difference approximations to our <5-derivatives. We see 
in Fig. H that (dP/dS) // = near T c for both values of 
ed- Dimerization increases the transfer integral (1 + 8) 



and tends to equalize charges on D and A sites then pro- 
ducing charge fluxes of opposite direction on the N and 
I side. (dP/dS) * ; = marks the point where charges are 
maximally delocalized and locates the N-I crossover in 
dimerized chains. 

Since Mott insulators on the I side are inherently corre- 
lated, all NIT features have been viewed in terms of cor- 
relations, including the divergence of dynamical charges 
at the NIT of the regular chain ]l0| . To address this 
point we investigated the non-interacting case, and found 
P(A, 6, U — 0) analytically for the infinite chain: 

P(A,^ = 0) = |--^§=ce% clP) a ! 6) (4) 

where cel{q Cl p, a, b) is the complete elliptic integral of the 
third kind M, calculated for q 2 c = (A 2 + AS 2 ) /(A 2 + 4), 
p = (A 2 + 4c5 2 )/(A 2 ), a = 1, and b = 0. For the regular 
chain the 5-derivative of p and of the elliptic integral 
above both vanish so that the frozen and the fluctuating 
charge contributions to 8Pj 88 are the coefficients of 5 in 
the first and second term of Eq. (||), respectively. For 
typical a values the frozen charge contribution smoothly 
increases from to ~ 0.05 as p increases from to 1. 
The fluctuating charge contribution is more interesting 
and is shown in the inset of Fig. || for positive A ( A < 
simply reverses the role of D and A) . (dP/ 85) /; diverges 
as A — * 0, much as it does for U — oo at NIT. Indeed, 
Eq. (g) at U = A = describes the half-filled band 
whose instability Peierls pointed out, and A = is the 
NIT for free electrons. We conclude that divergent IR 
intensity at the FE transition is not due to electronic 
correlations, but rather indicates a structural instability 
that displaces electronic charges. The behavior of chains 
with finite at the equilibrium dimcrization (dashed and 
dotted lines in the inset) in the N phase is similar to the 
U — oo case, with (dP/ 85) fi vanishing at A = and at 
r c , respectively. 

Fig. H shows results for Peierls coupling = 0.28 and 
0.64 in Eq. (|l|) with U = oo. The larger dimerization for 
6d = 0.64 strongly localizes electrons and reduces finite- 
size effects: rings of 10 and 12 sites appear superimposed, 
while they show small differences for =0.28. The IR 
intensity is dominated by fluctuating charges; it is small 
far in the N regime, peaks sharply at the FE transition, 
and decreases rapidly to vanish at the N-I crossover. In 
fact, since the frozen charge contribution is small, vanish- 
ing IR intensity can be taken as an operational definition 
of the charge crossover. The IR intensity smoothly in- 
creases again on the I side due to the increasing frozen 
charges and reaches typical values for ionic crystals as 
p — * 1. The peak in the IR intensity at the FE transi- 
tion is larger the nearer the transition is to the NIT of 
the regular chain, where electronic charges are most re- 
sponsive. IR intensity does not measure the strength of 
e-ph coupling, nor the dimerization amplitude, rather it 




FIG. 3: GS equilibrium ionicity, dimerization amplitude, 
polarization and IR intensity for U = oo and two different td- 
Note the scale change for both 8 and (dP/dS) 2 . Continuous 
and dashed lines refer to N=12 and 10, respectively (the two 
curves coincide in the left panels). 

is governed by the response of the electronic system and 
reflects its non-linearity. 

The Peierls-Hubbard model in Eq. ([!]) contains nei- 
ther intersite electron-electron interactions nor coupling 
to molecular (Holstein) modes. Consequently it cannot 
describe materials such as TTF-CA [jjj with a discontin- 
uous NIT. But Eq. (|l]) provides a useful starting point 
for systems with a continuous NIT, or, more precisely, 
a dimerization transition near a continuous NIT. Several 
CT salts dimerize in the N regime |o[ |§] with de- 
creasing temperature or increasing pressure; this mod- 
ulates r through the crystal's Madelung energy. Far 
IR data on these systems are however scanty, the only 
available results having recently been published on TTF- 
QB1-Q3 j22|. This compound has p ~ 0.3 at ambient 
conditions and dimerizes at T^ 68 K. Concomitantly, its 
ionicity increases from ~ 0.35 to 0.5 to reach ~ 0.6 at T 
~ 10K. Reflectivity measured at the lowest frequency (25 
cm -1 ) increases up to 90% at the FE transition, suggest- 
ing the proximity to a metallic state there. After dimcr- 
ization, the reflectivity decreases again. Approaching the 
FE transition, a progressive displacement of the IR inten- 
sity towards lower frequencies can be inferred from data 
in Ref. p^|, and the published IR intensity vs. T closely 
resembles the = 0.28 peak in Fig. ^[ Therefore we 
assign the amplified far IR intensity observed near the 
FE transition to the soft dimerization mode, and explain 
its abrupt drop after the transition as due to the van- 
ishing of the fluctuating contribution to the IR intensity 
at the N-I crossover. This interpretation differs from the 
one proposed in Ref. Q , where the far-IR absorption is 
assigned to the dynamics of neutral-ionic domain-walls. 
Additional experimental and theoretical work is needed 
to clarify this point. 

There are more data on static dielectric constants, £ 
pl[ p2| . Systems undergoing dimerization near a contin- 
uous NIT show a large peak in e just where long-range 
dimerization sets in. The similarity between the e peak 
(see e.g. Fig. 1 in ref. p2]|) at the FE transition and 
the IR intensity of the Peierls mode in Fig. 3 is sugges- 
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FIG. 4: The equilibrium dimerization and dielectric constant 
for U = and = 0.28. Dashed and continuous line refer to 
the purely electronic and to the total e, respectively. 



that, quite irrespective of electronic correlations, lin- 
ear e-ph coupling, harmonic phonons and delocalized 
electrons produce large and non-linear effects originat- 
ing from a complex interplay of electronic and nuclear 
degrees of freedom. 
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tive. E-ph coupling shifts oscillator strength from elec- 
tronic states to lattice vibrations, and this displacement 
to lower energy amplifies the system's response to a static 
field. 

To the first order in e^, the static response can be esti- 
mated as the sum of an electronic plus a vibrational con- 
tribution: e = 1 + (xei + Xvib)/so, with Xvib = tym/w °c 
ed{dP/dS) 2 , defined in terms of the vibrational transi- 
tion dipole moment and frequency p^ ]. Analytical ex- 
pressions of Xei are available for U — |24|, and Fig. 
4 compares the total (electronic+vibrational) e with the 
purely electronic term. The lattice parameters, required 
to fix the absolute value of e, are set to representative 
values for mixed stack crystals from Ref. 25 . Whereas 
the magnitude of the two terms is comparable, the sharp 
peak at the FE transition is due e-ph coupling. Again our 
interpretation differs from that in Ref. | p2[ , where the 
amplification of the dielectric constant at the FE transi- 
tion was related to the frequency of the far IR mode, via 
the Lyddane-Sachs- Teller relation. We instead underline 
the role of the dimerization mode in displacing electronic 
charges. Consistent with available data 0, |2^], we pre- 
dict large e peaks at FE transitions occurring at interme- 
diate p, i.e. for systems that, due to either small Peierls 
coupling or to disorder, undergo dimerization very close 
to the electronic (NIT) instability. Preliminary e cal- 
culations on U — oo open chains yield similar results, 
suggesting once more that the physics of FE transitions 
has more to do with phonons than with electronic corre- 
lations. 

Experimental evidence is accumulating for strong e-ph 
coupling in high-T c superconductors |jl|, || . The proxim- 
ity to an electronic instability is a unifying features of 
exotic superconductivity in cuprates, organics, or MgB2, 
as well as IR intensity or dielectric responses of CT salts. 
Using an adiabatic model for an extended electronic sys- 
tem with an electronic instability, we demonstrate here 
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